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Preface 


This  report  by  Dr.  John  Tweed  is  the  outgrowth  of  some  of  the  author's 
work  on  fracture  mechanics  initiated  while  he  was  a  Visiting  Assistant 
Professor  at  North  Carolina  State  University  in  1969-70. 

Although  the  major  portion  of  this  work  was  not  supported  by 
Grant  AFOSR-69-1779r  it  it»  being  presented  as  a  project  report  as  a  part  of 
our  continuing  interest  in  problems  of  fracture  mechanics. 

Dr.  Tweed  plans  to  obtain  some  computational  results  using  the  theoretical 
tools  derived  in  this  report. 


Project  Director 
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ABSTRACT 


In  this  paper,  the  author  makes  use  of  a  recent  development  in  the 
theory  of  Mellin  transforms  to  show  that  the  stress  intensity  factor  and  the 
crack  energy  of  a  crack,  which  originates  at  the  edge  of  a  circular  hole  in 
an  infinite  elastic  solid,  are  related  in  a  simple  fashion  to  the  solution  of 
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THE  DISTRIBUTION  OF  STRESS  NEAR  THE  TIP  OF  A  CRACK  WHICH 
ORIGINATES  AT  THE  EDGE  OF  A  CIRCULAR  HOLE 

by 

John  Tweed 


1.  INTRODUCTION, 

The  problem  of  determining  the  distribution  of  stress  near  the  tip  of 
a  crack  which  originates  at  the  edge  of  a  circular  hole  in  an  infinite  elastic 
solid  appears  to  have  been  considered  first  by  0.  L.  Bowie  [1]  who  solves  it 
by  using  a  complex  mapping  technique.  It  would  seem,  however,  that  the 
results  given  by  Bowie  are  not  very  accurate,  so  in  this  paper  we  wish  to 
show  that  the  stress  intensity  factor  and  crack  energy  are  related  to  the 
solution  of  a  Fredholm  equation  and  may  therefore  be  calculated  to  a  high 
degree  of  accuracy. 

We  shall  assume  that  the  problem  is  to  be  solved  under  the  conditions 
of  plane  strain  and  that  the  crack  and  the  hole  are  defined,  in  plane  polar 
coordinates  (r,  0),  by  the  relations  R  <_  r  £  Rb,  0=0  and 
0  r  <  R,  0  £  6  £  2tt  respectively. 


Figure  1. 
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If  the  loading  is  symmetric  about  the  plane  of  the  crack  it  is  clear 
that  the  problem  may  be  reduced  tc  that  of  finding  a  solution  of  the  equations 
of  elasticity  for  the  region  R  <  r  <  ®,  O<0  <  tt»  which  is  such  that 


(1) 

at  infinity  the  stresses  ^ri,(r* 

0),  ^(r,  6>r  cree^r'  0)  are 

0(r~2) 

and 

the  displacements  uy 

(r,  0),  u0(r,  0)  are  OCr”1), 

(2) 

°rS(r- 

0) 

*  0,  R  <  r  < 

(3) 

°r0(r' 

tt) 

•  u0(r,  u)  »  0,  R  <  r 

<  »» 

W 

■Vr<8' 

e) 

*  0,  0  <  0  <  Tt, 

(5) 

0re<R- 

e; 

*  0,  0  v  9  <  it, 

(6) 

°ee(r- 

0) 

*  -  f(r),  R  <  r  <  Rb, 

(7) 

us(r- 

0)  * 

0,  Rb  <  r  <  «, 

and 

(8) 

limit 

aue 

(T-°>  <  .  . 

r— *-R. 


dr 


2 -  REDUCTION  OF  THE  PROBLEM  TO  AN  INTECRAL  EQUATION. 

tn  order  to  find  a  suitable  representation  for  the  stresses  an!  dis¬ 
placements  in  the  problem  Get  out  above  we  shall  begin  by  superimporing  the 
solutions  of  problems  i  and  2  below. 

PROBLEM  1.  Find  a  solution  of  the  equations  of  elasticity  for  the  region 
R  r  <  <»,  0  <  e  <  ir9  which  is  such  that 

-2  -1 

(a)  at  infinity  the  stresses  are  0(r  )  and  the  displacements  are  0( r  ), 

(b)  ara(r,  0)  »  u0(r,  0)  «  0,  R  <  r  <  », 
and 

(c)  c  (r,  it)  3  u  (r,  k)  «  0,  R  <  r  <  *  . 

By  the  method  of  separation  of  variables  it  is  not  difficult  to  show 
(e.g,  see  Durelli,  Phillips  arid  T3ao  •°J)  that  the  Airy  stress  function 
for  this  problem  is  given  by 

^  or 

i(r,  C)  *  c  log  r  +  c,r  Accs  6  +  }  [c  r  ''  +  d  r  n^‘" 

o  1  %  n  n 

n=2 


jcoe  n0  , 


(2-1) 
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and  the  corresponding  stresses  and  displacements  by 

00 

o  (r,  e)  =  c  r  ^  -  y  [n(n  +  l)c  r  n  +  (n  +  2)(n  -  l)d  r  n]cos  n6  , 
rr  o  *• ,  n  n 

n=l 


-n-2 


cr  _ (r,  0)  =  -  l  [n(n  +  l)cnr“u  *  +  n(n  -  l)dQr  n]  sin  n0  , 
n=l  1  1 


and 


(2.2) 

(2.3) 


a._(r,  0)  =  -  c  r  2  +  l  [n(n  +  l)c  r  n  2  +  (n  -  2)  (n  -  l)d  r  n]cos  n0  ,  (2.4) 

00  o  _  ,  n  n 

n*i 


u  (r,  0)  *  / -  c  r  1  +  c,r  2cos  0 

r  F-  \  o  1 

00  « 

+  V  [ncrn^+(n+2-  4n)d  r  n+^]cos  n0\  , 
n-2  n  n  •  J 


(2.5) 


/  1  +  n  f  -2 

u0(r,  0)  =*  — g— 1  Jr  *  * 


00 

yc1r”2sin0  +  J  [n  c^r  n  1  +  (n  -  4  +  4ri)dnr  n+1]sin  n8>,  g) 
l  n=2  •* 


where  E  is  the  Young's  modulus  and  n  is  the  Poisson's  ratio  of  the  material: 
PROBLEM  2.  Find  a  solution  of  the  equations  of  elasticity,  for  the  half-plane 
0<^r<^®,  0  <  0  <  it,  which  is  such  that 

-2  -1 

(a)  at  infinity  the  stresses  are  0( r  )  and  the  displacements  are  0( r  ), 

(b)  at  the  origin  the  stresses  and  displacements  are  bounded, 

(c)  o  (r,  0)  =  0,  0  _<  r  <  », 
and 

(d)  arQ(r»  tt)  -  u0(r,  tt)  =  0,  0  _<  r  <  »  . 

By  utilising  the  properties  of  the  Mellin  transform  (e,g.  see  Tranter  [3]) 
it  can  be  shown  that  the  solution  of  the  problem  may  be  written  in  the  form 


o  (r,  0)  *  r  2M  "*■ 
rr 


A(s) 


2  stmts 

-  (s  +  2)cos( 


(s  +  4)cos(6  -  tt)  (s  +  2) 


(0  -  tr)s|  ; 


r 


(2.7) 
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/  \  -2,-1 
0re^r’  9)  =  r  M 


(s.  .+.  {sln(0  -  ir)(s  +  2)  -  sin(6  -  ir)s}  ;  r 

2  sln>  s 


i  2.8) 


°e6(r*  9) 


-2..-1 
r  M 


A(s) 


2  sinirs 


{ (8  +  2)cos(0  -  tt)s 


s  cos (0  —  tt)  (s  +  2)  >  ;  r 


(2.9) 


u  (r,  0)  * 


rE 


A(s) 


2(s  +  l)simrs 


{(s  +  2)cos(0  -  tt)8 


(s  +  4  -  4n)cos(8  -  it )  ( 8  +  2)}  ;  r 


(2.10) 


and 


u  (r,  6)  -  M'1 

6  rE 


Ms). 


2(s  +  l)8lnTT8 


{ (s  +  2)sin(0  -  t)s 


-  (s  -  2  +  4n)sin(8  -  tt)(s  +  2)}  ;  r 


-1 


(2.11) 


where  M  is  the  inverse  Mellin  transform  and  -1  <  Re(s)  <  0. 

Clearly  by  superimposing  the  solutions  of  these  two  problems  we  obtain 
a  solution  of  the  equations  of  elasticity  for  the  region  R  <  r  <  J  <  8  <  tt 
which  automatically  satisfies  conditions  (1),  (2),  and  (3)  and  which  is 
such  that 


o  (r,  0;  -  c  r  -  J  (n(n  +  l)c  r  +  (n  +  2)(n  -  l)a  r  ]cos  n0 
rr  o  n  n 

n=l 


,  -2..-1 
+  r  M 


Afs ) 


l(s  +  4)cos(6  -  tt)  (s  +  2) 


L2  sinns 


(s  +  2)cos(8  ■■  tt)s}  ;  r  , 


v<r' e)--  r, 

1  i 


-n-2 


n(s  t  r  +  n(n  -  l)d  r  sin  n0 
n  n  j 


,  -nl 
Jnr  1 


(2.1?; 


x  -2U-1 

+  r  M 


(s  t  2)A(s) 
2  suit's 


{sin(0  -  tt)(s  +  2)  -  sin(0  -  ir)s}  )  r 


(2. 13) 


5 


a  (r,  6)  =  -  c  r  2  +  l  [n(n  +  l)c  r  n  2  +  (n  -  2)(n  -  l)d  r  n]cos  n0 
oo  o  i  n  n 

n=l 


+  r“2*f 1  - 


A(s) 


2  sinus 


{(s  +  2)cos(6  -  ir)s 


-  S  C08 (0  -  "Xs 


+  2)}  ;  r 


v»r(r,  0) 


■L+  n 


-c  r'1  +  c.r_2cos0  +  V  {nc  r”n'*1  +  (n  +  2-~4n)d  r  n  1 
o  l  **„  n  n 

n»2 


rE 


MsL 


2(s  +  1)  sinus 


{ (s  +  2)cos(0  -  it) 8 


-  (s  +  4  -  4r})cos(0  -  ir)(s  +  2)}  ;  r 


and 


u0(r,  0)  - 


1  +  n 


-n-1 


-n+1, 


c1r~28ine  +  l  {nc  r_n  ^  +  (n  -  4  +  4n)dnr  n+1}sin  n0 
n»2  n  n 


+  xj_!ilir1 

rE 


.Alai 


2(8  +  1)  sinus 


{ (8  +  2)sin(0  -  tt)s 


-  (s  -  2  +  4ri)sin(0  -  u)(s  +  2)}  ;  r 


where  -1  <  Re(s)  •:  0  The  complete  solution  of  the  problem  may  now  be 
obtained  by  choosing  the  unknown  function  A(s)  and  the  unknown  sequences 
{c^}  and  {d^}  in  sucn  a  way  that  the  remaining  boundary  conditions  are 
satisfied. 

From  (2,14)  and  (2,16)  we  see  that  conditions  (6)  and  C7)  will  be 
satisfied  if  A(s)  is  a  solution  of  the  dual  equations 

M  ^(ACsHotus  ;  r)  =<  -  r2f(r)  -  r2F(r),  R  <  r  <  Rb 
M  *[(1  +  s)  ^AC's) ;r ]  *  0  ,  Rb  <  r  <  * 


(2.14) 

}cos  n0 

(2.15) 


(2.16) 


(2-17) 


6 


where  -1  <  Re(s)  <  0  and 


OO 

F(r)  *  -c  r  +  V  [n(n  +  l)c  r  n  +  (n  -  2)  (n  -  l)d  r  n]. 


o  t'1 

n=l 


If  we  now  assume  that  A(s)  may  be  written  in  the  form 


(2.18) 


rRb  . 

A(s)  -  p(t)  t8+1< 
*  R 


we  find  that  (see  Tweed  [4]) 


M”1[A(s) (1  +  s)"1;  r]  =  r  i 


0, 

Rb 

r 

0, 


dt 


0  <  r  <  R 
p(t)dt,  R  <  r  <  Rb 
Rb  <  r  <  ® 


and 


M  1[A(s)  cot-rrs  ;  r]  *  -  [  ^ 


Rb 


dt 


R  t  -  r 


and  hence  that  the  equations  (2.17)  will  be  satisfied  if 

-  [  L »  _  r  f(r)  _  r  F(r)  R  <  r  <  Rb. 
77  Jr  t  -  r 


(2.19) 


(2.20) 


(2.21) 


(2.22) 


The  equation  (2,22)  is  well  known  and  Tricomi  [5]  has  shown  that  its  solution 
is  given  by 


t  p(t)  - 


I  /  t  -  P 
Rb  -  t 


1/2 


r 

JR  \  y  -  R 


U/2 


y  f(y)  +  yF(y)  dy 

y  -  t 


(2.23) 


[(Rb  -  t)(t  -  R)] 


172  ' 


where  C  is  an  arbitrary  constant.  In  order  to  determine  C  we  make  use  of 

condition  (S'*  which  together  with  (2.16)  and  (2.20)  clearly  implies  that 

limit  p(r)  exis  a  and  hence  that  C  *  0.  It  now  follows  that 
r  — rR, 
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p(t)  is  given  by  the  expression 


1  I  t  -  R]1'2  j Rb  /Rb  -  y}1^2  y  f(y;  +  y  F(y) 

y  -  t. 


■  a  |5^ij  j 


R  \y  -  E 


dy. 


Similarly,  on  applying  conditions  (4>  and  (5)  we  find  that 


*  c  +  l  fn(n  +  1) R~ n  +  (n  +  2)(n  -  l)d  R~n+/i]cos  n9 


■*n+2. 


ri*l 


n 


=  M  1j — { (s  +  4) cos (0  -  ir)(s  +  2)  -  (s  +  2)cos(9  -  rr)s)  ;  R 
^2  sin-rrs 

00 

l  tn(n  +  l)c  R  n  +  n(n  -  l)d  R  n+2Jsin  n0 
n=l  n  n 


=  M 


-1 


^  ^  {sin(0  -  ir)(s  +  2)  -  sin(9  -  ir)s}  ; 


2  sitiTS 

where  0  <  e  <  rr,  and  hence  that 


c  -  -f-  M"1 lA(s> 
o  2tt 


-  2Rn  .-1 
c  3  -  M 

n  it 


{' 


s  +■  2  s 


s  s 


+  2/  *  R  ’ 

a(3)(s  +  2) 


(s2  -  n2) (s  +  2  +  n) 


;  R 


»  n  >,  1 


and 


d 

n  * 


_ A(s)  (s  +  2) 

♦  n)((s  +  2 ] 2  -  n2) 


5  R 


,  n  >  2 


where  -1  <  Re(s)  '  0< 

On  subs’’ iru ting  fr:m  (2:19)  into  (2,24)  through  (2  26)  and  working 
out  the  inverse  MelJm  transforms,  we  now  find  that 

-  R2  -Rb  -1 
■=  —  t  P(t)dt  » 


Rn  Rb  r  'p 

C  -S-  r  pi,,  }S.  -  j* 

r  2'.  JR  V>  ■  •  I- 


r.+2 


n  -  2 


I* 


n  [t/ 


dt ,  n  ■>  1 


(2,23) 


(2,24) 

(2-25) 

(2-26) 

(2-27) 


(2  28) 
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and 


,n-2  ,-Rb 


2it 


iR  ^  p«  {rM  i?r  -  i?n  *■  -i* 


(2.29) 


and  hence  that 


fRb 


F(r) 


-  f  p(t)  K(r,  t)  dt  , 

71  Jr 


where 


K<r,  t) . » V  -  it  .  iipuq  . 

t(R^  -  rt;  (R  -  rt)  R  -  rt 


J2  2 

t _  R  -  t 


2 

r  t 


1 

r 


(2.30) 


(2.31) 


At  this  point  we  find  it  convenient  to  introduce  the  function  P(t) 
which  in  defined  by  the  equation 


P(t)  =  [(t  -  R) (Rb  -  t)]1/2p(t). 


(2-32) 


On  substituting  from  (2.30)  into  (2.23)  and  taking  account  of  (2.32)  we  see 
that  P(t)  must  satisfy  the  integral  equation 


P(t)  - 


Rb 


P(p)  M(t,  p) 


R  [(p  -  R) (Rb  -  p)] 


1/2 


dp  *  S(t), 


where 


S(t)  = 


t  -  R 


’'t  J  R 


Rb 


Rfa  -  yl1/2  y  f(y)  dy 


y  -  R I 


y  -  t 


and 


M(tr  p)  = 


t  -  R 


it  t  ^  R 


■Rb 

f  Rb  -  y\ 

1 /2  . 

,  y  K(yt  p) 

Jr 

i  y  -  Rj 

y  -  t 

(2-33) 


(2.34) 


(2-35) 


If  we  now  substitute  from  (2-31)  into  (2.35)  we  find  that  M(t,  p)  may  be 
written  in  the  form 
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rr  t 


1/2 

(IT  -  pt)‘(p  -  R)] 


(t  -  R) (&  -  e)  Pt  (bp  -  R)J _ 

'  -2  .2,  „,l/2 


B(b  -  l)p 


2(R2  -  pt)(bp  -  R)1/2(p  *  R) 3/2 


1/2 


pt 


(t  -  R)R2(bp  -  R)1/2 
nt(R2  -  DtXe  -  R)1/2 


(2.38) 


3.  THE  STRESS  INTENSITY  FAt"  I  OR  AND  THE  CRACK  ENERGY. 

We  shall  now  show  that  the  stress  Intensity  factor  K  and  the  crack 
energy  W  which  are  defined  by  the  equations 

K  -  -  [2(Rb  -  t))1/2  ^  <r,0> 

r— *-Rb  2(1  -  n  ) 


(3.1) 


and 


W 


Rb 


o0(r,  0)  u0(r,  0)  dr 


(3.2) 


respectively,  are  simply  related  to  the  function  P(t)  which  was  introduced  in 
the  lest  section. 

By  substituting  from  (2.32)  into  (2.20)  and  taking  account  of  (2.16) 
we  see  that 


u0(r,  0) 


E 


Rb 


mLAL 


-  r 


[(t  -  R) (Rb  -  t)) 


1/2 


R  <  r  <  Rb. 


(3.3) 


It  now  follows  that  K  is  given  by 

yl/2 

K  - - fTT  p(Rb)  (3.4) 

[R(b  -  l)]lU 

and  W  by 


W 


E 


fRb 

Jr 


P(t)  dt _ 

((Rb  -  t)(t  -  R) ] 1/2 


•t 

f (r)  dr. 
R 


(3.5) 
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